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Most of the properties of the Earth-Moon system can be explained by a collision between a 
planetary embryo and the growing Earth late in the accretion process 1 3 . Simulations show 
that most of the material that eventually aggregates to form the Moon originates from the 
impactor 1 ®BD However, analysis of the terrestrial and lunar isotopic composition show them 
to be highly similai® n . In contrast, the compositions of other solar system bodies are sig¬ 
nificantly different than the Earth and Moon 12 14 . This poses a major challenge to the giant 
impact scenario since the Moon-forming impactor is then thought to also have differed in 
composition from the proto-Earth. Here we track the feeding zones of growing planets in a 
suite of simulations of planetary accretion 15 , in order to measure the composition of Moon¬ 
forming impactors. We find that different planets formed in the same simulation have distinct 
compositions, but the compositions of giant impactors are systematically more similar to the 
planets they impact. A significant fraction of planet-impactor pairs have virtually identical 
compositions. Thus, the similarity in composition between the Earth and Moon could be a 
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natural consequence of a late giant impact. 


Successful models for Moon formation typically require a relatively low-velocity, oblique 
impact 1 between the proto-Earth and up to a fewxO.l Earth-mass {Me) planetary embryo. Such 
Moon-forming impacts typically occur at the late stages of planetesimal accretion by the terrestrial 
planet^ 3 1 A circum-terrestrial debris disk is formed from material ejected during these impacts. 
The composition of the disk, in a typical impact, is dominated by material from the impactor 
mantle (> 60 weight percent, % wf 1 4 5 ) with a smaller contribution (typically « 20%) from the 
proto-Earth. More material can be extracted from the proto-Earth when a slightly sub-Mars sized 
body hits a fast spinning planet, that is later slowed down by resonances 16 . The spin should be 
close to the break-up velocity. Another possible channel producing a significant mixed material 
from both the planet and impactor is the rare collision between two comparable mass embryos, 
which both masses are of about half of Earth’s mass 17 . Although these new models can potentially 
solve some of the composition issues borne by the giant-impact scenario, they do require ad hoc 
assumptions and pose several difficulties (see ref. 3 for a discussion). Here, we focus on the former 
typical giant-impact events, in which the Moon aggregates mostly from material originating from 
the impactor. 

Lunar meteorites and rock samples returned by the Apollo mission have a very similar com¬ 
position to Earth’s mantle, across a variety of different isotopesP 11 . Combining these with the 
giant-impact simulation results, one infers that the Moon-forming impactor and the Earth should 
have had a similar composition. This poses a fundamental difficulty to the giant-impact model for 
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the origin of the Moon, since analysis of material from other solar system bodies have shown them 
to significantly differ from that of the Earth (see refs. 12 13 and ref. 14 for a review). This would 
suggest that the composition of the Moon-forming impactor should have similarly differed from 
that of the Earth, in contrast with the giant-impact basic prediction. 

Here, we analyze the results of extensive N-body simulations of terrestrial planet formation 
to show that the Earth-Moon composition challenge can be alleviated. In particular, we show that 
the compositions of a significant fraction, 20% to 40%, of giant-impactors are consistent with being 
similar to that of the planets they impact. More generally, late giant impactors have significantly 
more similar compositions to the planets they impact compared with other planets in the same 
system, showing large differences. 

To study the compositions of planets and their impactors we analyzed 40 dynamical simula¬ 
tions (from ref. 15 ; using the Mercury code®) of the late stages of planetary accretion, following the 
formation of Jupiter and Saturn, and after all the gas in the protoplanetary disk has been dissipated 
and/or accreted to the gas-giants (see ref. 19 for a recent review). Each simulation started from 
a disk of 85-90 planetary embryos and 1000-2000 planetesimals extending from 0.5 to 4.5 AU. 
Jupiter and Saturn are fully formed and have different orbits and inclinations in different sets of 
simulations (detailed descriptions can be found in ref. and in the Methods. Within 100-200 mil¬ 
lion years, each simulation typically produced 3-4 rocky planets formed from collisions between 
embryos and planetesimals. Each of these planets accreted a large number of planetesimals during 
its evolution. All collisions were recorded and provide a map of each planet’s feeding zone (see 
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also ref. 20 ). Assuming that the initial composition of material in the protoplanetary disk is a func¬ 
tion of its position in the initial protoplanetary disk, one can compare the compositions of different 
bodies formed and evolved in the simulations. 

Previous studies explored the compositions of the different planets formed in similar sim¬ 
ulations. However, the composition of impactors on formed planets have been hardly explored. 
Pahlevan & Stevenson 21 analyzed a single statistically limited simulation which included a total of 
~150 particles, and compared the compositions of any impactors on any planets during the simula¬ 
tion (not only giant impacts, due to small number statistics) to that of the planets. They concluded 
that the scatter among the compositions of the various impactors is comparable to the observed 
differences between the planets. In particular, they found that none of the planetary impactors in 
the simulation they analyzed had an isotopic composition similar enough to the final planet to yield 
an Earth-Moon-like composition similarity. 

Using the data from our large set of high resolution simulations, we compare the composition 
of each surviving planet with that of its last giant impactor, i.e. the last planetary embryo that im¬ 
pacted the planet (typical impactor-to-planet mass ratio in the range « 0.2 — 0.5; see Table 1). We 
only include the 20 cases where both the impactor and planet are composed of at least 50 particles 
each, so as to have significant statistics. Analysis of the additional data for impactors composed 
from a smaller number of particles (and hence smaller statistics for the specific composition) are 
consistent with the higher resolution cases discussed here, as shown in the Methods. The compar¬ 
ison is then done as follows. First we compare the feeding zones of the planet and impactor, as 
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shown in the examples in Fig. 1 (the cumulative plots for these and all other cases can be found 
in the Methods). We calculate the probability (P) that the feeding zones of the impactor and the 
planet are drawn from the same distribution, using a two-group Kolmogorov-Smimov (KS) test 
(probabilities shown in the plots and in Table 1). In 3 out 20 cases the feeding zones contributing 
to the Moon and those contributing to the planet are consistent with being drawn from the same 
parent distribution. In other words, the Moon feeding zones, if derived solely from the impactor, 
are consistent with the Earth’s in 15% of the impacts. The consistency further improves if we as¬ 
sume that a fraction of the proto-Earth was mixed into the Moon (as suggested by detailed collision 
simulations showing a 10 — 40% contribution from the proto-Eartrf 14 ). For the typical 20% mix 
of proto-Earth material with the impactor material forming the Moon (as found in simulations), 
35% of cases are consistent, and the success rate increases further for a higher mass contribution 
from the proto-Earth (see Table 1 and Extended Data Figures 1 and 2). While this shows that the 
proto-Earth and the Moon-forming impactor may have had similar feeding zones, it does not yet 
quantitatively guarantee that the composition is as similar as that of the Earth-Moon system. 

We therefore further explore the compositional similarities with the Earth-Moon system, and 
calculate the oxygen isotope ratios of our simulated planets. We assume that a linear gradient 
existed in the 1 ' O isotopic composition in the initial protoplanetary disk of the solar-system. Fol¬ 
lowing Pahlevan & Stevenson 21 we calibrate the initial 17 0 isotopic composition in each of our 
simulations using Earth and Mars’ measured compositions (see ref. 21 and the Methods for more 
details; where we also discuss the sensitivity of results to the calibration used, as well as the cri¬ 
teria for which planet-impactor pairs are considered in the analysis. We find qualitatively similar 
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results when using different criteria and calibrations, as we discuss in detail in the Methods). Given 
this calibration we assign each planetesimal and planetary embryo a specific initial 17 0 isotopic 
abundance based on its initial orbit, and then average the contribution of all accreted planetesimals, 
while weighting each accreted planetesimal/embryo according to its mass, to obtain the 17 0 of the 
planets and impactors and derive the offset between them. 

The 17 0 isotope is chosen for comparison since it provides the most stringent constraint 
on the Earth-Moon similarity, and its abundances were measured across a variety of solar system 
bodies (enabling the best opportunity for calibration). The measured difference between Mars’ and 
Earth’s 17 0 abundances (used for calibration) is A 17 OMars = +321 ± 13 parts per million (ppm) 11 
(a similarly large difference was found for the composition of 4 Vesta asteroid derived from HED 
meteorites; —250 ± 60 ppm). The difference between the Earth and Moon is just A 1 ' O = 12 ± 3 
ppm! 11 ! In Table 1 we show the A 11 0 differences between the impactors and the planets in our 
simulations. Note that we adopt the same Earth-Mars composition difference calibration as used by 
Pahlevan & Stevenson 21 . The results are linearly dependent on the adapted calibration; see Table 
1. The calibration factor is defined by C ca i = (A 17 0 4 — A 1 '0 3 )/A 17 0 M ars 5 where A 17 0 4 and 
A 1 '0 3 refer to the compositions of the fourth and third planet in the simulations (unless only three 
planets formed, in which case the third and second planet were taken for calibration purposes). 

In 20% of the cases the impactors and planets have absolute offsets comparable or smaller 
than the measured absolute offset in the Earth-Moon system, i.e. smaller than the lcr limit es¬ 
timated using the Lunar samples (< 15ppm). Taking into account the la uncertainty calculated 
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for the A 17 O in the simulated systems, the fraction of consistent pairs can be as large as 40% of 
the whole sample. This fraction becomes larger when partial mixing of Earth material is allowed, 
as observed in simulation data (it increases to 50% (55%), for 20% (40%) contribution from the 
planet; see Table 1 and Fig. 2). Even planet-impactor pairs with statistically different feeding 
zones have A 1 ' O offsets significantly smaller than those found for Mars and Vesta, in most cases. 
More generally, planet-impactor pairs are robustly more similar in composition compared with 
pairs of surviving planets in the same system (see Fig. 2, Extended Data Figures 3, 4 and 5, as 
well as the Supplementary Information Table, for the A 1 ' O difference distribution for planets and 
impactors). Not less important, the differences between the planets are of similar order to those 
found between the Earth, Mars and Vesta, i.e. consistent with the observations of the solar system. 
Interestingly, a small fraction of the planets do have very similar composition (small A 1 ' O differ¬ 
ence), suggesting the possibility for the existence of solar system bodies with similar compositions 
to the Earth besides the Moon. As shown in Extended Data Table 2 and in the Methods this result 
still holds when considering lower thresholds for the minimal number of particles composing the 
planet/impactor (between 1 to 40). In particular, the mean fraction of compatible planet-impactor 
pairs extends between 10% and 20% for all cases. The fractions become even higher (20%-40%) 
when accounting for the la uncertainties. 

The Earth-Moon composition similarity poses a major challenge to the standard model of the 
giant-impact scenario as it conflicts with the predominant derivation of the Moon composition from 
the mantle of the impacting planet 14 . It therefore gave rise to a wide range of alternative impact 
scenarios 14 E3 17 21 24 . However, all of these models suffer some potentially considerable difficulties 


7 




and/or require fine-tuned conditions (see ref. 14 for a review). Our analysis of solar-system-like 
planet formation scenarios potentially offers a solution to the major composition-similarity ob¬ 
stacle, for the standard giant-impact scenario. We find that a significant fraction of all planetary 
impactors could have had similar composition to the planets they impacted, in contrast with the 
composition of different planets existing at the same planetary system. Note that the solution of 
impactor-planet similar composition suggested by our results well applies for the origin of the 
A 1 2 '0 similarity between the Earth and the Moon, and may similarly apply for the other isotopic 
similarities of the Silicon and Tungsten 25 . However, it is still debated whether even a similar 
impactor-planet composition could resolve the composition similarity of the latter (e.g. Silicon), 
as Earth’s silicate mantle may reflect the consequences of silicon sequestration by a core formed at 
high temperatures on a large planetary body 25 2 ^, i.e. larger than the typical impactors considered. 

We conclude that our findings can potentially resolve the apparent contrast between the ob¬ 
served similarity of the Earth and the Moon composition and its difference from that of other solar 
system bodies. This primordial composition similarity solution may therefore lift the prime ob¬ 
struction for the standard giant-impact origin of the Moon, as well as ease some of the difficulties 
for alternative giant impact scenarios suggested in recent years. 
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Model 

# 

M P 

(M 0 ) 

Mi 

(M 0 ) 

M l 

Adp 

N P 

Ni 

^coll 

(Myrs) 

C C alA 17 0 (ppm) 

KS-probability 

0% 

20% 

40% 

0% 

20% 

40% 

cjsl5 

1 

0.94 

0.43 

0.46 

123 

97 

50.7 

13 ± 14 

10 ± 13 

8 ± 12 

0.0039 

0.13 

0.67 

cjsl5 

2 

0.78 

0.27 

0.35 

209 

78 

80.9 

(-1.05 ± 0.26) X 10 2 

-84 ± 23 

-63 ± 22 

1.3 X 10 -7 

3.1 X 10 -4 

0.079 

cjsl 

3 

1.25 

0.42 

0.34 

219 

73 

149.7 

64 ± 13 

52 ± 11 

39.9 ± 9.6 

5.3 X 10 -11 

8.1 X 10 -7 

0.0092 

cjsl 

4 

1.05 

0.39 

0.37 

128 

78 

186.2 

(-1.97 ± 0.20) X 10 2 

(-1.57 ± 0.20) X 10 2 

(-1.18 ± 0.20) X 10 2 

1.1 X 10 -29 

5.1 X 10 -18 

0.052 

cjsl 

D 

1.21 

0.38 

0.31 

219 

75 

123.5 

-24 ± 17 

-20 ± 15 

-15 ± 13 

0.0023 

0.056 

0.19 

cjsecc 

6 

0.94 

0.36 

0.38 

117 

79 

80.4 

-51 ± 34 

-40 ± 30 

-31 ± 27 

3.4 x 10~ 4 

0.025 

0.33 

cjsecc 

7 

1.01 

0.32 

0.32 

144 

68 

75.5 

-12 ± 16 

-10 ± 13 

-7 ± 12 

0.038 

0.050 

0.20 

cjsecc 

8 

1.02 

0.42 

0.41 

148 

89 

36.9 

13 ± 79 

11 ± 71 

8 ± 66 

0.054 

0.21 

0.53 

eejsl5 

D 

0.70 

0.19 

0.27 

111 

52 

77.1 

9.1 ± 7.4 

7.3 ± 6.2 

5.5 ± 5.2 

0.071 

0.18 

0.32 

eejsl5 

10 

0.55 

0.13 

0.24 

263 

65 

24.6 

98 ± 29 

78 ± 24 

59 ± 22 

1.6 X 10 -12 

2.3 X 10 -8 

0.0025 

eejsl5 

11 

0.78 

0.22 

0.29 

256 

69 

102.3 

(-1.08 ± 0.35) X 10 2 

-87 ± 31 

-65 ± 27 

2.9 X 10 -13 

7.9 X 10 -7 

0.0069 

eejsl5 

12 

0.73 

0.26 

0.36 

298 

87 

105.6 

26 ± 18 

21 ± 16 

16 ± 14 

1.3 X 10 _1 ° 

5.5 X 10 -6 

0.034 

eejsl5 

13* 

1.30 

0.33 

0.25 

525 

126 

199.8 

93 ± 1.5 X 10 2 

74 ± 1.3 X 10 2 

55 ± 1.1 X 10 2 

7.1 X 10 -8 

5.4 x 10 -6 

0.0043 

eejsl5 

14 

0.50 

0.18 

0.36 

170 

53 

33.1 

-26 ± 75 

-21 ± 65 

-16 ± 56 

0.14 

0.13 

0.24 

eejsl5 

15 

0.50 

0.13 

0.26 

234 

61 

32.0 

-73 ± 56 

-58 ± 48 

-44 ± 43 

2.7 X 10~ 8 

1.8 X 10 -7 

0.062 

eejsl5 

16 

0.67 

0.33 

0.49 

213 

120 

145.0 

(1.37 ± 0.51) X 10 2 

(1.10 ± 0.45) X 10 2 

83 ± 42 

2.4 x 10 -5 

0.013 

0.43 

eejsl5 

17* 

1.15 

0.41 

0.36 

177 

69 

168.3 

(7.9 ± 1.3) X 10 2 

(6.3 ± 1.1) X 10 2 

(4.75 ± 0.93) X 10 2 

3.5 X 10 -19 

1.2 x 10 _1 ° 

0.0028 

ejsl5 

18 

0.81 

0.36 

0.44 

63 

55 

76.3 

-81 ± 48 

-65 ± 42 

-49 ± 37 

1.5 X 10 -5 

0.0039 

0.29 

jsres 

19 

1.04 

0.32 

0.31 

166 

89 

79.5 

55 ± 22 

44 ± 19 

33 ± 17 

3.0 X 10 _1 ° 

1.5 X 10 -6 

0.0078 

jsres 

20 

1.27 

0.63 

0.50 

134 

89 

176.8 

-84 ± 29 

-67 ± 25 

-50 ± 23 

5.6 X 10 -7 

0.011 

0.50 

Obs. 




A 17 O (ppm) 


Earth 


1 


0 ± 3 


Moon 


0.012 


12 ± 3 


79 Mars 


0.07 


321 ± 13 


4 Vesta 


4.33 X 10~ 5 


-250 ± 80 



* 3-planet systems; calibration was done on the 1 st and 2 nd planets. 


Table 1: Properties of the modeled planet-impactor systems, and comparison with the observations of Solar system bodies. Mp,Np 

and MpNi arc the mass and number of particles in the planet and the impactor, respectively; t co \\ is the collision time in the simulations; C ca iis 
the calibration pre-factor (see main text). The A 17 0 composition difference and the KS-probability (for the planet and impactor feeding-zone 
distribution to be sampled from the same parent distribution) arc shown both the case of contribution of planetary material to the newly-formed 
Moon, and the cases of 20% and 40% contribution of material from the planet. 
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Figure 1: The distribution of planetesimals composing the planet and the impactor. Panel a 
shows a case where the origins of the planetesimals composing the planet and the impactor are 
consistent with being sampled from the same parent distribution for the expected typical 20% 
contribution of planetary material in moon-forming impacts (KS-test probability > 0.05). Panel b 
shows a case where the planet and impactor compositions are inconsistent (P < 0.05), but become 
consistent once a significant (40%) contribution of material from the planet is considered. The 
lower plots in each panel show the results assuming a different contributions from the planet (four 
cases are shown 10%, 20%, 30% and 40%). The cumulative distribution for these cases as well as 
all other planet-impactor pairs in Table 1 can be found in the Methods. 
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Figure 2: The cumulative distribution of the absolute A 17 0 differences between planets and 
their last giant impactors (blue), compared with the differences between planets in the same 
system (red). Panels a, b and c correspond to the cases of zero, 20% and 40% contribution of 
material from the planet to a Moon formed from these impacts, respectively. The vertical lines 
depict the A 1 ' O difference of the Earth-Moon system (dashed lines for the ±cr). The differences 
between the planet-impactor pairs are systematically smaller than those found between different 
planets (the same parent distribution for the two groups can be excluded with high confidence; KS 
probability 6.7 x 10~ 8 ,1.1 x 10 -8 and 1.3 x 10 -9 for the zero, 20% and 40% cases, respectively). 
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Methods 


In the following we supply additional data on the planet-impactor pairs’ composition and the com¬ 
positions of different planets at the same systems. We also provide a more detailed information on 
the methods used as well as discuss the sensitivity of our results to the various criteria and calibra¬ 
tions which we applied. The full simulations data used in this work can be provided by the authors 
upon request. 

Initial configuration of Jupiter and Saturn 

The simulations analyzed here are described in detail in ref. ^ 15 . These various simulations explore 
a range of different initial conditions for the gaseous planets. In particular, in the cjs and cjsecc 
simulations Jupiter and Saturn are placed on orbits with semimajor axes of 5.45 and 8.18 AU and 
mutual inclination of 0.5°. In cjs the orbits are circular while in cjsecc they are eccentric with 
ej = 0.02 and e$ = 0.03. In eejs Jupiter and Saturn are placed on their current position (5.25 
and 9.54 AU), with mutual inclination of 1.5° and larger eccentricities than the observed ones 
(ej = es = 0.1 or ej = 0.07 and es = 0.08). In ejs the orbits of Jupiter and Saturn have similar 
parameters to those observed (aj = 5.25AU and ej = 0.05, as = 9.54 and es = 0.06) with 
mutual inclination of 1.5°. Finally, in jsres, Jupiter and Saturn are placed at aj = 5.43AU and 
as = 7.30AU with ej = 0.07 and es = 0.01 and with mutual inclination of 0.2°. 
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The composition difference between planets in the same system 


The Supplementary Information Table shows the A 17 0 differences between the different planets 
(in each system) and the impacted-planets analyzed in the main text. The full cumulative distri¬ 
bution for these data can be seen in Fig. 2 in the main text. Note that in those cases where two 
impacted planets were analyzed in the same system, the differences are shown both in respect to 
the first and second planets. 

The cumulative composition distribution of planet-impactor pairs 

We used the following procedure to calculate the spatial distribution of the feeding zones (used 
for Fig. 1, Extended Data Fig. 1 and 2 and Table 1 KS-probabilities). We extracted the record of 
planetesimals that constitute the planet and the impactor before the last moon-forming impact, as 
well as the planet composition after the collision. In order to account for the different contribution 
coming from particles of different masses we replicated n* times each particles, where n, is the 
ratio between the mass of the i-Th. particle and the minimum mass of the planetesimals. The plan- 
etesimal record was then used to produce the distribution of the feeding zones used in our analysis. 
In cases where contribution of planetary material to the Moon composition was considered, we 
randomly chose particles from the planet and added them to the planetesimals composing the im¬ 
pactor, where appropriate numbers of particles were taken so as to produce the relevant fractional 
contribution (for the different cases of 10, 20, 30 or 40% contribution). We then repeated the same 
analysis as done for the impactor, with these new mixed impactors. 
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A 17 O calibration 


In order to calculate the A 1 '0(=5 1 '0 — 0.52<5 16 O) for the planet and impactor pairs we followed 
the procedure described by ref. 21 . In order to assign specific values of A 17 0 to each particle in the 
simulation we calibrated our simulations with the solar system observations. We assume a linear 
gradient of A 17 0 with heliocentric distance r 


A 17 0(r) = cir + c 2 . (1) 

where the two free parameters in Equation 1 were calibrated imposing that the third planet formed 
in the system has the composition of the Earth (A 1 ' O = 0%o) and the fourth one has the compo¬ 
sition of Mars (A 17 0 = +0.32%o). In cases where only three planets formed (marked with a * in 
Table 1), we assigned the composition of the Earth and Mars to the second and third planets in the 
simulation, respectively. We then mass-averaged over all the A 17 0(r) of the planetesimals that 
accreted forming the Earth. We used the initial position of each body as the heliocentric distance 
r. We did the same for the planetesimals composing “Mars”. In this way we have the system of 
equations 


J2 m i.E(c'iri.E+ c 2) 

Me 

N ' 

J2iiH.m(civ uM +C2} 

Mm 


= 0 

= 0.32, 
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where r i)E and m, E (r l M and m r hI ) are the initial position and mass of the i th planetesimal com¬ 
posing the Earth (Mars) and M E (M m ) is the final total mass of the Earth (Mars). In this way it 
has been possible to evaluate the A 1 ' O value for each planetesimal in the system and thus for all 
the planets in each system. Given this calibration we evaluate the A 1 '© of the planet (A 17 Op) 
and of the last impactor (A 17 0/), as the average of the A 1 ' O of all their respective components, 
as well as calculated the la SEM for each of these values (o> and a{). To check whether or not 
the planet-Moon system is consistent with the Earth-Moon system we evaluated the difference 
A 1 '0/ — A 1 ' Op and the relative error a = \J a 2 P + erf. 

In order to calculate the A 17 0 when a fractional contribution from the planet is included, we 
added the average A 1 ' O of the planet and the impactor, each weighted according to the appropriate 
fractional contribution considered. 

In order to study the sensitivity of our results to the Earth-Mars composition difference cal¬ 
ibration used, we also considered lower and higher calibrations, between 0.5-1.5 times the Earth- 
Mars A^O difference. We re-analyzed the fraction of compatible planet-impactor pairs (i.e. pro¬ 
ducing planet-moon pairs with composition difference equal or smaller than the Earth-Moon com¬ 
position difference) for these different calibrations. The results are summarized in Extended Data 
Table 1. We find that although the fraction of consistent pairs decrease with the use of larger dif¬ 
ference calibration, as expected, difference calibrations as much as 1.5 larger than the Earth-Mars 
difference still give rise to a mean 5% of planet-impactor pairs with similar composition (and 40% 
within la uncertainty in the simulation compositions), rising to a mean 10% to 20% for the cases 
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of 20% and 40% mixing of the planetary material, respectively. In other words, the results are 
generally robust to that level and do not dependent on a fine-tuned calibration. 

Dependence on the criteria for the planet-impactor pairs considered in the analysis 

In order to verify the robustness of our results to different criteria for the choice of planet-impactor 
pairs used in our analysis we studied various different criteria: (1) Use all planet - last-impactor 
pairs, considering smaller thresholds (i.e. smaller than the 50 particles threshold considered in the 
main text) for the number of composing particles (and corresponding masses), but requiring an 
impactor mass of at least 0.5 M Mars to assure a moon-forming impact. Taking a threshold of 1, 
10, 20, 30 and 40 minimal number of particles, we find that the general conclusion is unchanged; 
impactors have more similar composition to the planets they impact compared with other plan¬ 
ets in the system. The mean fraction of planet-impactor pairs with comparable similarity as the 
Earth-Moon system is between « 10 — 20% in all cases (and up to 20 — 40% considering the 
1 a uncertainties), as shown in Extended Data Table 2. Extended Data Figures 3, 4 and 5 show 
the cumulative distributions of the composition of planets and last impactors for all the systems, 
regardless of the number of particles that contributed to their formation, and for minimum of 10, 
20, 40 and 50 particles composing the planet and last impactor. The Extended Data Figures 3, 4 
and 5 are shown for 0%, 20% and 40% mixing between planet’s and impactor’s material. 

(2) Consider only last-impactors on the third planet. Once we require the impactor to have at least 
0.5 M Mars and be composed of a significant number of particles the statistics become too small. 
The only case for which we have a significant (18 cases) is when we do not consider a minimal 
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threshold for the number of composing particles. In this case we find 2 out of 18 (11%; up to 
pa 30 % when with the 1 a uncertainty in the simulations composition) planet-impactor pairs with 
composition difference equal or smaller than the Earth-Moon system. 
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Factor 

0% 

20% 

40% 

Mean 

la 

Mean 

la 

Mean 

la 

0.5 

35% 

50% 

35% 

60% 

35% 

70% 

0.75 

20% 

50% 

25% 

50% 

35% 

60% 

1 

20% 

40% 

20% 

50% 

25% 

55% 

1.25 

5% 

40% 

20% 

50% 

20% 

50% 

1.5 

5% 

40% 

10% 

40% 

20% 

50% 


Extended Data Table 1: The mean fraction of last impactors with compatible composition 
respect to the planet they impact is given for different normalization factors and mixing 
percentages. The pairs which are consistent within la of the simulation uncertainties are also 
given. 
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N ■ 

0% 

20% 

40% 

N 

1 v cases 

Mean 

la 

Mean 

la 

Mean 

la 

0 

10.1% 

21.3% 

10.1% 

25.8% 

11.2% 

33.7% 

89 

10 

10.5% 

24.6% 

10.5% 

29.8% 

12.3% 

36.9% 

57 

20 

12.5% 

30% 

12.5% 

35% 

15% 

42.5% 

40 

30 

14.3% 

35.7% 

14.3% 

42.9% 

17.9% 

46.4% 

28 

40 

18.2% 

36.4% 

18.2% 

45.5% 

22.7% 

50% 

22 

50 

20% 

40% 

20% 

50% 

25% 

55% 

20 


Extended Data Table 2: The mean fraction of planet-impactor consistent pairs is shown for 
different mixing percentages (0, 20, 40%) and minimum numbers (N min ) of particles com¬ 
posing the impactor and planet (1,10, 20, 30, 40, 50). The pairs which are consistent within la 
of the simulation uncertainties are also given. Only the N cases cases in which the last-impactor has 
a mass > 0.5M Ma rs have been taken into account. 
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Extended Data Figure 1: The cumulative distribution of the planetesimals composing the 

planet and the impactor, showing all planet-impactor pairs in Table 1, cases 1-12 (panels a- 

1), including the cumulative distributions corresponding to the histograms in Fig. 1 in the main 
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text. 


































































Extended Data Figure 2: The cumulative distribution of the planetesimals composing the 
planet and the impactor, showing planet-impactor pairs in Table 1, cases 13-20 (panels m-t). 
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Extended Data Figure 3: The cumulative distribution of the absolute A 17 0 differences between 
planets and their last giant impactors (blue), compared with the differences between planets 
in the same system (red) assuming 0% mixing between Earth and Moon material. From 
the top left panel (a) to the bottom right panel (f) we consider all the systems, regardless of the 
number of particles that contributed to their formation, and planets and last impactors composed 
by a minimum of 10, 20, 40 and 50 particles. Only last-impactors with mass > 0.5M Mars have 
been taken into account. 
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Extended Data Figure 4: The cumulative distribution of the absolute A 17 0 differences between 
planets and their last giant impactors (blue), compared with the differences between planets 
in the same system (red) assuming 20% mixing between Earth and Moon material. From 
the top left panel (a) to the bottom right panel (f) we consider all the systems, regardless of the 
number of particles that contributed to their formation, and planets and last impactors composed 
by a minimum of 10, 20, 40 and 50 particles. Only last-impactors with mass > 0.5M Ma rs have 
been taken into account. 
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Extended Data Figure 5: The cumulative distribution of the absolute A 17 0 differences between 
planets and their last giant impactors (blue), compared with the differences between planets 
in the same system (red) assuming 40% mixing between Earth and Moon material. From 
the top left panel (a) to the bottom right panel (f) we consider all the systems, regardless of the 
number of particles that contributed to their formation, and planets and last impactors composed 
by a minimum of 10, 20, 40 and 50 particles. Only last-impactors with mass > 0.5MMars have 
been taken into account. 
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Supplementary Information 


Model 

Diff. # 

Mass 

a( AU) 

Ref. Planet 

A 17 0 (ppm) 

cjsl5 

1 

3.02 x HT 1 

0.498 

cjsl5 #1.1 

-27 ±21 

cjsl5 

2 

9.40 x 10' 1 

0.742 

cjsl5 #1.2 

0± 13 

cjsl5 

3 

1.45 

1.41 

cjsl5 #1.3 

(-2.82 ±0.21) x 10 2 

cjsl5 

4 

6.10 x 10- 2 

2.06 

cjsl5 #1.4 

(-6.02 ±0.79) x 10 2 

ejsl5 

5 

3.68 x 10" 2 

2.35 

cjsl5 #1.5 

(-5.85 ±0.49) x 10 2 

cjsl5 

6 

7.83 x 10" 1 

0.580 

cjs 15 #2.1 

0 ± 25 

cjsl5 

7 

5.57 x 10" 1 

0.891 

cjs 15 #2.2 

-77 ± 30 

cjsl5 

8 

9.78 x HT 1 

1.37 

cjs 15 #2.3 

-287 ±31 

cjs 15 

9 

5.76 x HT 1 

1.91 

cjs 15 #2.4 

-606 ± 46 

cjsl 

10 

1.25 

0.720 

cjsl #3.1 

0 ± 7.2 

cjsl 

11 

1.05 

1.94 

cjsl #3.2 

-231 ± 19 

cjsl 

12 

8.80 x 10- 2 

2.33 

cjsl #3.3 

-109 ±31 

cjsl 

13 

2.47 x 10" 2 

2.45 

cjsl #3.4 

-429 ± 23 

cjsl 

14 

1.25 

0.720 

cjsl #4.1 

231 ± 19 

cjsl 

15 

1.05 

1.94 

cjsl #4.2 

0 ± 25 

cjsl 

16 

8.80 x 10" 2 

2.33 

cjsl #4.3 

122 ± 35 

cjsl 

17 

2.47 x 10" 2 

2.45 

cjsl #4.4 

-198 ± 28 

cjsl 

18 

1.21 

0.765 

cjsl #5.1 

0± 14 


29 



























































































































cjsl 

19 

1.32 

1.83 

cjsl #5.2 

406 ± 28 

cjsl 

20 

2.25 x 10- 2 

3.55 

cjsl #5.3 

86.4 ±9.3 

cjsecc 

21 

6.85 x 10” 1 

0.574 

cjsecc #6.1 

128 ± 32 

cjsecc 

22 

9.35 x 10- 1 

0.901 

cjsecc #6.2 

0 ± 28 

cjsecc 

23 

5.83 x 10" 1 

1.41 

cjsecc #6.3 

-211 ±33 

cjsecc 

24 

1.61 x HT 1 

2.45 

cjsecc #6.4 

-531 ± 36 

cjsecc 

25 

1.01 

0.557 

cjsecc #7.1 

0± 12 

cjsecc 

26 

7.18 x HT 1 

0.953 

cjsecc #7.2 

72 ± 12 

cjsecc 

27 

5.93 x HT 1 

1.52 

cjsecc #7.3 

-49 ± 14 

cjsecc 

28 

4.99 x 10- 2 

1.98 

cjsecc #7.4 

(-3.69 ±0.43) x 10 2 

cjsecc 

29 

4.86 x 10- 2 

2.36 

cjsecc #7.5 

96 ±39 

cjsecc 

30 

3.21 x HT 2 

2.51 

cjsecc #7.6 

-320.9 ±8.3 

cjsecc 

31 

2.21 x HT 1 

2.76 

cjsecc #7.7 

-283 ± 52 

cjsecc 

32 

1.02 

0.622 

cjsecc #8.1 

0 ± 38 

cjsecc 

33 

5.61 x HT 1 

1.13 

cjsecc #8.2 

-228 ± 46 

cjsecc 

34 

6.92 x HT 1 

1.38 

cjsecc #8.3 

-585 ± 47 

cjsecc 

35 

1.40 x HT 1 

2.12 

cjsecc #8.4 

-905 ± 52 

cjsecc 

36 

3.59 x 10" 2 

2.81 

cjsecc #8.5 

(1.0 ±1.2) x 10 2 

eejsl5 

37 

6.97 x HT 1 

0.565 

eejsl5 #9.1 

0 ± 4.5 

eejsl5 

38 

8.15 x 10” 1 

0.905 

eejsl5 #9.2 

-70.8 ±6.7 
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eej s 15 

39 

1.58 x 10" 1 

1.61 

eejsl5 #9.3 

-79.9 ±9.4 

eejsl5 

40 

5.20 x 10- 2 

3.15 

eejsl5 #9.4 

-398.9 ±3.2 

eejsl5 

41 

5.47 x 10” 1 

0.540 

eejsl5 #10.1 

0 ± 20 

eejsl5 

42 

7.77 x lO' 1 

0.852 

eejsl5 #10.2 

(3.01 ±0.24) x 10 2 

eejsl5 

43 

1.06 x HT 1 

1.20 

eejsl5 #10.3 

(8.69 ±0.24) x 10 2 

eejsl5 

44 

3.44 x 10” 1 

1.51 

eejsl5 #10.4 

(5.48 ±0.47) x 10 2 

eejsl5 

45 

5.00 x 10" 2 

1.61 

eejsl5 #10.5 

(8.41 ±0.23) x 10 2 

eejsl5 

46 

5.47 x HT 1 

0.540 

eejsl5 #11.1 

(-3.01 ±0.24) x 10 2 

eejsl5 

47 

7.77 x 10" 1 

0.852 

eejsl5 #11.2 

0 ± 27 

eejsl5 

48 

1.06 x 10" 1 

1.20 

eejsl5 #11.3 

(5.67 ±0.78) x 10 2 

eejsl5 

49 

3.44 x 10" 1 

1.51 

eejsl5 #11.4 

(2.47 ±0.49) x 10 2 

eejsl5 

50 

5.00 x 10" 2 

1.61 

eejsl5 #11.5 

(5.40 ±0.26) x 10 2 

eejsl5 

51 

7.34 x HT 1 

0.590 

eejsl5 #12.1 

0± 14 

eejsl5 

52 

6.77 x 10" 1 

0.917 

eejsl5 #12.2 

-35 ± 14 

eejsl5 

53 

4.39 x HT 1 

1.64 

eejsl5 #12.3 

-110 ± 16 

eejsl5 

54 

5.92 x 10" 2 

3.10 

eejsl5 #12.4 

-430 ±17 

eejsl5 

55 

6.06 x 10" 2 

3.92 

eejsl5 #12.5 

-604.8 ±9.7 

eejsl5 

56 

1.30 

0.654 

eejsl5 #13.1 

0 ± 99 

eejsl5 

57 

2.97 x HT 1 

1.37 

eejsl5 #13.2 

(1.6 ±1.4) x 10 2 

eejsl5 

58 

3.16 x 10” 1 

1.75 

eejsl5 #13.3 

(-1.6 ± 1.5) x 10 2 
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eej s 15 

59 

3.34 x 10" 1 

0.518 

eejsl5 #14.1 

(2.25 ±0.70) x 10 2 

eejsl5 

60 

5.68 x 10” 1 

0.692 

eejsl5 #14.2 

164 ± 56 

eejsl5 

61 

4.95 x 10" 1 

0.925 

eejsl5 #14.3 

0 ± 54 

eejsl5 

62 

2.11 x 10” 1 

1.26 

eejsl5 #14.4 

(-3.2 ±0.70) x 10 2 

eejsl5 

63 

2.26 x 10" 1 

1.59 

eejsl5 #14.5 

(-1.37 ±0.80) x 10 2 

eejsl5 

64 

3.70 x 10' 1 

2.62 

eejsl5 #14.6 

(-1.360 ±0.038) x 10 3 

eejsl5 

65 

4.97 x HT 1 

0.538 

eejsl5 #15.1 

0 ± 42 

eejsl5 

66 

4.64 x 10" 1 

0.748 

eejsl5 #15.2 

(1.40 ±0.41) x 10 2 

eejsl5 

67 

6.72 x HT 1 

1.08 

eejsl5 #15.3 

(8.03 ±0.44) x 10 2 

eejsl5 

68 

9.02 x 10- 2 

1.62 

eejsl5 #15.4 

(4.83 ±0.63) x 10 2 

eejsl5 

69 

8.73 x 10" 2 

2.01 

eejsl5 #15.5 

(1.39 ±0.10) x 10 3 

eejsl5 

70 

4.97 x HT 1 

0.538 

eejsl5 #16.1 

(-8.03 ±0.43) x 10 2 

eejsl5 

71 

4.64 x 10" 1 

0.748 

eejsl5 #16.2 

(-6.63 ±0.42) x 10 2 

eejsl5 

72 

6.72 x HT 1 

1.08 

eejsl5 #16.3 

0 ± 45 

eejsl5 

73 

9.02 x 10" 2 

1.62 

eejsl5 #16.4 

(-3.20 ±0.64) x 10 2 

eejsl5 

74 

8.73 x 10" 2 

2.01 

eejsl5 #16.5 

(5.9 ±1.0) x 10 2 

eejsl5 

75 

1.15 

0.646 

eejsl5 #17.1 

0 ± 85 

eejsl5 

76 

6.74 x 10" 1 

1.27 

eejsl5 #17.2 

(-1.96 ±0.91) x 10 2 

eejsl5 

77 

7.31 x 10- 2 

1.67 

eejsl5 #17.3 

(-5.2 ±2.7) x 10 2 

ejsl5 

78 

4.73 x HT 1 

0.506 

ejsl5 #18.1 

(3.04 ±0.32) x 10 2 
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ejsl5 

79 

7.93 x 10" 1 

0.808 

ejsl5 #18.2 

(2.39 ±0.29) x 10 2 

ejsl5 

80 

8.08 x HT 1 

1.38 

ejsl5 #18.3 

0 ± 35 

ejsl5 

81 

5.31 x lO" 2 

2.37 

ejsl5 #18.4 

(-3.20 ±0.50) x 10 2 

ejsl5 

82 

1.27 x 10” 1 

3.04 

ejsl5 #18.5 

(-8.07 ±0.44) x 10 2 

jsres 

83 

1.04 

0.621 

jsres #19.1 

0± 14 

jsres 

84 

1.27 

1.12 

jsres #19.2 

113 ±17 

jsres 

85 

1.01 x HT 1 

2.17 

jsres #19.3 

(5.67 ±0.62) x 10 2 

jsres 

86 

3.42 x lO" 2 

2.41 

jsres #19.4 

(2.47 ±0.33) x 10 2 

jsres 

87 

1.48 x HT 1 

2.44 

jsres #19.5 

(3.07 ±0.77) x 10 2 

jsres 

88 

7.92 x lO" 2 

2.57 

jsres #19.6 

(7.59 ±0.52) x 10 2 

jsres 

89 

1.57 x 10” 1 

2.7672 

jsres #19.7 

(1.74 ±0.20) x 10 2 

jsres 

90 

1.70 x 10" 1 

2.7654 

jsres #19.8 

(5.92 ±0.62) x 10 2 

jsres 

91 

1.04 

0.621 

jsres #20.1 

-113 ± 17 

jsres 

92 

1.27 

1.12 

jsres #20.2 

0 ± 20 

jsres 

93 

1.01 x HT 1 

2.17 

jsres #20.3 

(4.52 ±0.63) x 10 2 

jsres 

94 

3.42 x lO" 2 

2.41 

jsres #20.4 

(1.34 ±0.34) x 10 2 

jsres 

95 

1.48 x HT 1 

2.44 

jsres #20.5 

(1.94 ±0.78) x 10 2 

jsres 

96 

7.92 x 10" 2 

2.57 

jsres #20.6 

(6.45 ±0.52) x 10 2 

jsres 

97 

1.57 x 10” 1 

2.7672 

jsres #20.7 

61 ±21 

jsres 

98 

1.70 x 10” 1 

2.7654 

jsres #20.8 

(4.78 ±0.62) x 10 2 
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Supplementary Information Table: The A 1 ' O differences between the planets in each analyzed 
system and the impacted-planets. See the main text and Methods for further details. 
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